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ABSTRACT 

We have used the Owens Valley Radio Observatory Millimeter Array to obtain 
a map of the J = (1^0) CO emission from the circumstellar shell around HD 
179821, a highly evolved G-type star which will probably explode as a supernova 
in the next 10^ yr. Very approximately, the gas presents as a circular ring with an 
inner diameter of 3'.'95, an outer diameter of ~12" and with azimuthal variations 
in the CO brightness by about a factor of 2. Until about 1600 years ago, the 
star was a red hypergiant losing about 3 x 10~^ Mq yr"^ at an average outflow 
speed of 32 km s~^. We propose that when HD 179821 explodes as a supernova, 
it may resemble Kepler's supernova remnant and thus some of the anisotropies in 
supernova remnants may be intrinsic. If the factors which cause the anisotropic 
mass loss in HD 179821 persist to the moment when the star explodes as a 
supernova, the newly-born pulsar may receive a momentum "kick" leading to a 
space motion near ~700 km s^^. Independent of the angular asymmetries, the 
radially detached shell around HD 179821 may be representative of environments 
which produce dust echoes from gamma-ray bursts. 

Subject headings: circumstellar matter -stars: mass loss - stars: winds-outflows 
- stars: supernovae -ISM: supernova remnants 



1. INTRODUCTION 



After main sequence 0-type stars of more than 20 Mq with T^jf > 30,000 K consume 
their interior hydrogen, they evolve to become red hypergiants with mass loss rates near 
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10-^ M0 yr-i (see, for example, Chin & Stothers 1990, Schaller et al. 1992). Ultimately, the 
star explodes as a type II supernova. Because deviations from spherical symmetry might be 
crucial in supernova explosions (Burrows 2000), we hope to gain insights into the physics of 
supernovae by studying the anisotropic circumstellar shells of pre-supernovae. 

The morphology of a supernova remnant is determined both by the history of mass 
loss from the pre-supernova and the structure of the supernova explosion itself. We hope 
to learn whether nonspherical shapes in supernova remnants result from nonspherical mass 
loss by the pre-supernova star, from intrinsically nonspherical supernova explosions, from 
irregularities in the surrounding interstellar medium, or from some combination of these 
different effects. 

A mystery related to supernova explosions is that the newly-born pulsars sometimes 
have space motions over 700 km s^^ (Cordes & Chernoff 1998). It is possible that these 
motions result from a "momentum-kick" acquired during a non-spherical explosion, another 
reason for studying asymmetries in the outflows from pre-supernovae. 

Gamma-ray bursts may be produced during supernova explosions. If so, then dust at 0.1 
to 1 pc from the star may result in a "dust echo" some 20-30 days after the initial burst (Esin 
& Blandford 2000). The morphology of dust around massive stars may provide constraints 
on models for dust echoes. 

The winds from red hypergiants typically contain large amounts of dust, and therefore 
these objects arc strong infrared sources (sec, for example, Jura & Klcinmann 1990). As 
discussed both by Jura & Werner (1999) and below, it is likely, though controversial, that HD 
179821 (my = 7.9 mag, b = -4.96°, spectral type G5Ia, and a strong IRAS source with F^(25 
/im) — 650 Jy) is a highly evolved star near the Humphreys-Davidson limit with an initial 
mass near 30 Mq that will become a supernova within the next ~10^ yr. The alternative 
model is that it was a ^^1 solar mass star that has become a pre-Planetary Nebula, and it lies 
at a distance of about 1 kpc from the Sun. Although dust nebula are observed around post 
main sequence hypergiants with luminosities greater than 10^ Lq (de Jager 1998, Nota et 
al. 1995, Voors et al. 2000), most of these stars lie in the Galactic Plane where interactions 
with the local interstellar medium are important. Because HD 179821 is removed from the 
Galactic Plane, it may be a circumstellar shell with relatively little contamination by swept- 
up interstellar gas. About ~20% of 0-type stars are runaways (Blaauw 1993) with speeds 
greater than 50 km s~^, and these stars can travel as much as 1000 pc from the Galactic 
Plane before they explode. We propose that HD 179821 is post-main sequence, runaway, 
massive star. 

Previously, we have reported Keck observations of the 11.7 iim emission from HD 179821 
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and found that the circumstellar shell is very approximately a circular ring with an inner 
diameter of about 3'/5 (Jura & Werner 1999). However, both the 11.7 fim emission and the 
map of the OH maser spots in the circumstellar envelope of HD 179821 (Claussen 1993) 
display significant deviations from circular symmetry. A major difficulty with interpreting 
the 11.7 fim map is that at this wavelength, hu/k = 1230 K, while the grains are typically 
colder than 140 K. As a result, the emission is from the Wien portion of the Planck curve 
and small uncertainties in the temperature translate into large uncertainties in estimating 
the mass. As with most bright evolved stellar infrared sources, HD 179821 also is a strong 
CO source (Zuckerman & Dyck 1986, Bujarrabal, Alcolea & Planesas 1992, van der Veen 
et al. 1993, Knapp et al. 2000, Josselin Sz Lebre 2001). Therefore, in order to achieve a 
better tracer of the mass loss history and kinematics of the outflow, we have acquired a high 
resolution map of the circumstellar CO emission with the Owens Valley Radio Observatory 
Millimeter Array. 

In order to learn more about the evolution of supernovae and their remnants, we compare 
our maps of the shell around HD 179821 with the appearance of Kepler's supernova remnant 
which, with an age of 400 years, is one of the youngest known remnants in the Milky Way 
(Burrows 2000). Since it lies at b = 6.8° and is at a distance from the Sun of ~ 4800 pc 
(Reynoso & Goss 1999), it is ~500 pc from the Galactic Plane where there is relatively 
little interstellar matter. Kepler's remnant thus presents an excellent opportunity to study 
the interaction of a supernova explosion with its pre-supernova circumstellar shell. Kepler's 
supernova remnant probably (Decourchelle & Ballet 1994, Rothenflug et al. 1994, Hughes 
1999) though not certainly (Kinugasa & Tsunemi 1999) resulted from the explosion of a 
massive star. 

In Section 2, we summarize our picture of HD 179821 and present an additional argument 
for why we think that it is a massive star. In Section 3, we present our observations. In 
section 4, we present a simple model to explain the data while in Section 5 we compare the 
asymmetries in the circumstellar matter around HD 179821 with those in Kepler's supernova 
remnant. In Section 6, we note how our data may help both understand why pulsars may 
have space motions near 700 km s~^ and the dust echoes from gamma-ray bursts. We present 
our conclusions in Section 7. 



2. THE DISTANCE AND LUMINOSITY OF HD 179821 

The distance to HD 179821 is not well determined by trigonometric parallax (the value 
measured with the Hipparcos satellite was 0.18 ± 1.12 x 10~^ "), and indirect means must be 
used. Two models have been proposed. (1) The star might be at a "large" distance of about 
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6000 pc, have a luminosity near 3 x 10^ L© and thus have had an initial main sequence mass 
near 30 Mq (Hawkins et al. 1995). (2) The star might lie at a "small" distance of about 
1 kpc, have a luminosity near 10^ Lq and be the post- Asymptotic Giant Branch (AGB), 
pre-Planetary Nebula descendant of a star that was initially ~1 Mq on the main sequence 
(Hrivnak, Kwok & Volk 1989). While inferences of the luminosity of the star from the 
atmospheric abundances might resolve this argument; to date they are controversial (Reddy 
& Hrivnak 1999, Lobel & Dupree 2000, Thevenin, Parthasarathy & Jasniewicz 2000), and 
not evidently conclusive. The detection of CO and Na I emission at 2.2 //m (Hrivnak, Kwok 
& Geballe 1994, Oudmaijer et al. 1995) suggests that standard model atmospheres may not 
apply, and additional methods to estimate the distance should be investigated. 

Jossehn & Lebre (2001) report a well cahbrated measure of the ratio of integrated 
intensities in the J = (2— >^1) lines of ^^CO and ^^CO of 3.2 from which they argue that HD 
179821 is a post-AGB star. However, Bujarrabal et al. (1992) reported a value of 9.6 for 

this ratio, and 6.6 for the ratio of the integrated intensities of the ^^CO and ^^CO J = (1^0) 
lines. Given that the hypergiant IRC+10420 has a value of ^^CO/^^CO of 9 ± 2 (Fix & Cobb 
1987) and that the CO rotational lines are optically thick, observations of the circumstellar 
CO isotope ratio do not provide a compelling method for determining the luminosity and 
origin of the star. 

From the properties of the circumstellar matter, we argue that HD 179821 is a massive, 
distant star. Zuckerman & Dyck (1986) first noted that the high LSR velocity of the star 
(+100 km s~^) and high average outfiow velocity of the gas (32 km s~^) are best understood 
if the star lies at ~6 kpc. There is a rough correlation between the outflow velocity, \4o, 
from a red giant and its luminosity (Jones, Hyland & Gatley 1983). Therefore, since Voo for 
HD 179821 is near 32 km s~^ (see below), which is significantly larger than the typical value 
for a mass-losing red giants of 15 km s~^, it is plausible that this star has a relatively high 
luminosity. Barnbaum, Kastner & Zuckerman (1991) found in a sample of 124 AGB carbon 
stars that 4 stars have wind outfiow speeds greater than 30 km s~^. These 4 objects lie at 
low galactic latitudes and presumably have main sequence progenitor masses over 2.5 M0. 
If HD 179821 lies at a distance of 1 kpc, then its high LSR radial velocity indicates that it 
is a member of an old population with a main sequence progenitor mass less than 1 Mq, an 
inference not consistent with its relatively high outfiow velocity. 

In addition to comparing \4o of HD 179821 with mass- losing red giants, we can also 
compare its outflow velocity with the measured values for Planetary Nebula. Approximately 
10% of all Planetary Nebulae have expansion speeds greater than 30 km s^^ (Sabbadin 1984, 
Weinberger 1989). However, such speeds appear to be a consequence of acceleration of the 
wind material after the star evolves beyond being a mass-losing red giant, since Gesicki & 
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Zijlstra (2000) find that expansion speeds greater than 30 km s~^ are usually found only for 
Planetary Nebulae with radii larger than 0.03 pc. For a presumed distance from the Sun of 1 
kpc, the current radius of the peak density around HD 179821 is 0.008 pc. We conclude that 
if HD 179821 is a pre-Planetary Nebula with an outflow speed greater than 30 km s~^ and 
the densest material lying within 0.008 pc of the star, then it is very unusual and perhaps 
unique. 

The relatively high value of the LSR velocity for HD 179821 is more easily understood 
if it is a distant object and its speed is a result of Galactic rotation. Figure 1 shows a plot of 
LSR velocity vs. CO outflow velocity for a complete sample of high mass-loss rate (> 10~^ 
Mq yr~^) oxygen-rich AGB stars that are within 1 kpc of the Sun and further north than 
S = -32° (Jura & Kleinmann 1989). We also show on this plot the location of HD 179821 
{I = 35.6, b = -5.0) and the well known hypergiant IRC +10420 {I = 47.0, b = -2.5, D = 
5 kpc, L = 5 X 10^ Lq, Jones et al. 1993) which hes relatively close to HD 179821 in the 
sky. The better agreement of the kinematics of the outflow from HD 179821 with those from 
1RC+ 10420, rather than with the outflow kinematics of the local AGB stars suggests that 
HD 179821 is a red hypergiant. If HD 179821 is a massive, runaway star, then its radial 
velocity cannot serve as an accurate measure of its distance. However, if it does lie at 6 kpc 
from the Sun, then at its Galactic longitude, it lies beneath the 5 kpc ring, the region in 
the disk of the Milky Way with the greatest concentration of interstellar CO and giant H II 
regions (Scoville & Sanders 1987). Therefore, if HD 179821 is a hypergiant, a distance of 6 
kpc from the Sun is plausible but not certain. 

Another argument for the high luminosity of HD 179821 is that it has a near-infrared 
reflection nebulae similar to that of the hypergiant IRC-l-10420 (Kastner & Weintraub 1995). 
Furthermore, Jura & Werner (1999) point out that the OH maser map presented by Claussen 
(1993) is more naturally understood if the star hes at a distance of 6000 pc rather than 1000 
pc. That is, if HD 179821 lies at "only" 1000 pc from the Sun, then a uniquely high ultraviolet 
circumstellar opacity is required to account for the spatial extent of the OH maser emission 
which is determined by the depth of penetration of the ambient interstellar radiation field 
into the outer regions of the circumstellar shell. 

Finally, based on its kinematics and historic data, we present another argument that HD 
179821 probably lies further than 1 kpc from the Sun and thus probably has a distance near 
6000 pc. The average expansion speed of 32 km s^^ of the circumstellar matter is determined 
from the profile of the circumstellar CO emission and is therefore independent of distance. 
In contrast, the inferred physical diameter of the inner boundary of the circumstellar is 
derived from both the measured inner angular diameter of 3'.'5 and the assumed distance. If 
HD 179821 hes at 6000 pc from the Sun, the time required for the gas in the inner shell to 
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expand to its current size is 1600 yr. Therefore, there would not have been much change in 
the star during the past 100 yr. On the other hand, if the star has just left the AGB and 
lies at a distance of 1000 pc, then the time to expand to its current size is only 260 yr, and 
during the past 100 yr, it could have undergone observable changes. Such historic changes 
of the light from post-AGB stars in fact have been measured. For example, during the past 
100 years, the Egg Nebula, a famous post-AGB star, has brightened from mB= 15 mag to 
m^ = 13 mag (Gottlieb & Liller 1976). However, there is no historic evidence for any major 
changes in either the brightness or spectral type of HD 179821. In the Henry Draper catalog, 
which is nearly 100 years old, the spectral type is G5, the same spectral type given in the 
SIMBAD data base from modern observations. Since the change in spectral type from G5 
either to G8 or to G2 corresponds to a change in the effective temperature of the star of 
~250 K (Drilling & Landolt 2000), then it appears that during the last 100 years, dT^ff/dt 
< 3 K yr~^ for HD 179821. In order to have evolved from the AGB when it was losing a large 
amount of mass with an effective temperature of 2600 K to its current effective temperature 
of either 6750 K (Reddy & Hrivnak 1999) or 5660 K (Thevenin et al. 2000) within 260 years, 
then the contraction of the star and consequent rise in its surface temperature should have 
proceeded with an average value of dT^ff/dt of >11 K yr~^, which is larger than inferred. 
Furthermore, the magnitude of 8.1 given in the Bonner Durchmusterung catalog, which was 
compiled around 1860, is approximately the same as its current average visual magnitude of 
7.9; there is no evidence for major changes during the past 140 years. Therefore, the historic 
evidence supports the view that HD 179821 is a distant, luminous star. 



3. OBSERVATIONS OF THE CIRCUMSTELLAR CO 

We used the six-element Owens Valley Radio Observatory MiUimeter Array during 
November - December 1999 to observe the J = (1^0) ^^CO emission from HD 179821. Phase 
calibration was performed on the quasar 1749+096 while flux calibrations were derived from 
3C 273 and 3C 454.3. The synthesized beam was 2'.'0 x 1'.'35 at position angle -85°; the la 
rms noise was 40 mJy/beam. The map of the integrated CO intensity is shown in Figure 
2 where the contour levels are spaced by 1.8cr with each level corresponding to 1.07 K or 
32.5 mJy beam~^. The channel maps with 5.35 km s~^ resolution are shown in Figure 3. 
The channel map velocities are given relative to the assumed LSR velocity of 100 km s"^ 
(Bujarrabal et al. 1992). The hne profiles are not symmetric, and estimates of hne center 
have ranged 95 and 105 km s~^ (Reddy & Hrivnak 1999; Josselin & Lebre 2001). The position 
of the star was taken from the Hipparcos data as q;(2000.0) = 19'* 13"* 58.61^ 5(2000.0) = 
00° 07' 31'.'93 and is shown on our maps. 
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We show in Figure 4 the azimuthally- averaged intensity as a function of angular offset 
from the star in the channel at line center and at a velocity offset of +27 km s~^. As can 
be seen from Figure 4, the peak of the CO emission is seen to have an angular diameter 
of 3'.'95, slightly larger than the 3'/5 diameter found in the 11.7 fim. image (Jura & Werner 
1999). In our CO data, the peak to the north is about a factor of 1.4 times brighter than the 
peak to the south; in the 11.7 yum image, the contrast between the northern and southern 
peaks is about a factor of 2 (Jura & Werner 1999). In the east- west direction, the shell is 
less well defined. As shown in Figure 4, the CO emission can be traced to a radius of about 
6" so that the diameter of the CO emission is at least 12". The total flux measured with 
the interferometer in this 12" region is 45% - 55% of the flux measured with the 30m IRAM 
telescope with a beam diameter of 22" (Bujarrabal et al. 1992, Jossclin & Lebre 2001). 
Therefore, the CO is somewhat more extended than shown in our data. Figure 5 shows the 
CO spectrum at the center of the circumstellar envelope. 



4. CIRCUMSTELLAR GAS DISTRIBUTION AROUND HD 179821 

4.1. A first approximation: spherical geometry 

In order to estimate the mass loss rate when the star was a red giant, we start with 
the simplification of spherical symmetry. An estimate of the gas loss rate during the red 
hypergiant phase, M (Mq yr~^), can be derived from the single-beam measurements of the 
circumstellar CO emission. According to Kastner (1992), for observations of the J = (1—^0) 
transition of CO with the 30m IRAM telescope, if M < 10~^ Mq yr^^ then 

M = 6xl0-'T^,{0fV^D^ (1) 

where Tmb{0) is the main beam antenna temperature (K) at line center, \4o is the gas outflow 
velocity (km s^^), and D is the distance (kpc) of the star from the Sun. In this expression, we 
assume that in the circumstellar envelope, [C0]/[H2] = 10^^, consistent with the abundance 
analysis of the photosphere of HD 179821 by Rcddy & Hrivnak (1999) and the assumption 
that when the star was losing mass as a red hypergiant, as much carbon and oxygen were 
incorporated into CO as possible. With T^b = 1.4 K (the average of the values reported by 
Josselin & Lebre 2001 and Bujarrabal et al. 1992), V^o — 32 km s~^ and D — 6 kpc, then 
M — 3 X 10~^ Mq yr~^. This inferred value of M is so large that the models by Kastner 
(1992) are not particularly accurate, but, as noted below, our CO maps are consistent with 
this result. For a steady-state wind with V^o = 32 km s~^ and M = 3 x 10"'^ Mq yr""'', 
n(i?2) = 1-4 X 10^ cm~^ where Rn denotes the distance from the star (10^^ cm). 

Jossehn & Lebre (2001) report recent single-beam observations of the CO emission in 
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the J = (1^0) and J = (2— >1) transitions. Their measured integrated intensities agree 
within 20% with those found by Bujarrabal et aL (1992) for the -^^CO emission, but they 
disagree by a factor of 3 for the ^^CO hne. Using essentially the same integrated intensity for 
the ^^CO line as we do, and assuming a distance of 1 kpc, Josselin & Lebre (2001) estimate a 
mass loss rate of 2 x 10~^ M© yr~^. Scaled to a distance of 1 kpc, we would estimate a mass 
loss rate of 8 x 10~^ Mq yr~^, a factor of 25 smaller than that found by Josselin & Lebre 
(2001). There are three contributions to this difference. First, Jossehn & Lebre (2001) use 
a formula very similar to equation (1) but with a numerical coefficient that is a factor of 2 
larger to estimate the mass loss rate from the CO (J = (1^0) line profile. Their equivalent 
to equation (1) is derived by Loup et al. (1993) based on a scahng procedure rather than 
on detailed theoretical calculations as employed by Kastner (1992). Second, we assume that 
[C0]/[H2] = 10-3 while Josselin & Lebre (1993) adopt [C0]/[H2] = 5 x 10-^ another factor 
of 2 in the estimate of the mass loss rate. Third, Josselin & Lebre (1993) implicitly assume 
that the CO emission is spatially limited and does not fill the telescope beam. This leads 
to the requirment that there must be more gas in a smaller volume and thus a factor of 6 
increase in their estimate of M over our inferred value. In view of the theoretical calculations 
by Mamom et al. (1988) and our map, we see no need to make this adjustment. 

With M > 10~^ Mq yr"^, the gas density is sufficiently high that the rotational levels 
of the CO molecule are mainly populated by collisions with H2 (Kastner 1992). Thus, we 
can interpret our maps of the CO emission by comparing them with models for coUisional 
excitation of CO which have been developed for interpreting observations of interstellar 
matter (Castets et al. 1990, Sakamoto 1996). In the outer envelope around HD 179821 
where n{H2) < 10^ cm-^, the CO is subthermally excited so that the population in the J = 
1 rotational level is lower than predicted by the Boltzmann weighting function. In this case, 
a collisional excitation of CO into the J = 1 level is followed by radiative de-excitation so the 
observed Rayleigh- Jeans brightness temperature, Tg, of the J = (1^0) CO emission depends 
directly upon the collisional excitation rate. Consequently, Tg scales as n(H2)^(7(T^) where 
g{TK) is a slow function of the kinetic temperature, Tk- For the likely case that TV > 20 K 
(see below), g{TK) varies as T^^^ (Castets et al. 1990). 

We now consider the azimuthally-averaged CO intensity at 4" offset (or Rn = 3.6) from 
the star, a region where the density is probably low enough {n^H^) ~ 1100 cm~^ from above) 
that this simple model for collisional excitation pertains, yet also a location where our data 
have a relatively high signal to noise. At this offset of 4", Tb(0) the brightness temperature 
in the channel at km s~^, is 5.5 K or 160 mJy/beam (see Figure 4). To reproduce this 
observation, we assume [C0]/[H2] = 10"^, Tk = 60 K (see below) and a spherically sym- 
metric outflow. In this case, for a sightline whose impact parameter is Rimpact, dV/dR, the 
velocity gradient used in the calculations by Castets et al. (1990), is approximately given by 
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Voo/ Rimpact- With these parameters, n(H2) = 600 cm ^ a,t R — Rimpact (Castets et al. 1990). 
This estimated density at Rii = 3.6 is in reasonable agreement with the prediction from the 
model derived above with M = 3 x 10"^ Mq yr'^ and V^o = 32 km s"^ that n(H2) = 1100 
cm~^. Therefore, within a factor of 2, it seems that an average mass loss rate of 3 x 10~^ 
Mq yr~^ can explain both the observed CO integrated line profile (Bujarrabal et al. 1992) 
and the maps presented here. 

Additional evidence in favor of the model for subthermal collisional excitation of the 
CO is given by the angular variation of Tg. In the outer envelope of the nebula, the model 
predicts that Tb(0) scales as n{Ri,npactY Rimpact- If M and \4o are constant, then n{Rimpact) 
varies as Ri^pacv thus Tb{0) is predicted to vary as where (j) is the offset angle from 
the central star. As shown in Figure 4, in the channel at km s""^, the CO intensity falls 
from 270 mJy/beam (0 = 3") to 40 mJy/beam = 6"). This decrease of Tb as is in 
reasonable agreement with the expectation of the model for subthermal collisional excitation. 
Furthermore, if Tb varies as then between an inner radius, (pin and and outer radius, (pout, 
the total flux varies as {(pout - <Pin)- With 0j„ = 2", the approximate measured value of the 
inner radius of the CO ring, and (pout = 6" and 11" for the OVRO interferometer and IRAM 
30m telescope measurements, respectively, then it is expected that F(OVRO)/F(IRAM) = 
4/9 or 0.44 in agreement with the measured values of 0.45 to 0.55 reported above. 

With M = 3 X 10"'' M0 yr~\ then MVoo equaled 6 x 10^^ g cm s"^ when the star 
was losing a large amount of mass. Currently, L^,/c equals 4 x 10^*^ g cm s"^. This ap- 
proximate agreement between MV^o and L*/c is consistent with models for winds driven by 
radiation pressure on dust (Lamers & Cassinelli 1999), observations of other mass-losing red 
hypergiants (Jura & Kleinmann 1990) and expectations for the pre-supernova evolution of 
massive stars (Heger et al. 1997). Also, this mass loss rate derived from the gas of 3 x 10~^ 
Mq yr~^ compares with the value of 4 x 10~^ Mq yr~^ derived by Jura & Werner (1999) 
from the emission by the circumstellar dust with an assumed gas to dust ratio of 100. 

The observed diameter of the CO emission of more than 12", which is substantially 
larger than the diameter of the OH emission of ~ 4-5" (Claussen 1993), can be understood 
in the usual models for the photo-dissociation of molecules by ambient interstellar ultraviolet 
radiation as they flow out of the star (Glassgold 1996). The OH molecule is mainly protected 
by dust and therefore can be photodissociated at ~2" from the star (Jura & Werner 1999). 
The CO is self-shielding, and for a model with M = 3 x 10"'^ M© yr"^ and Voo = 30 km s~^, 
the predicted fractional abundance is reduced as a result of photodissociation by a factor 
of 0.5 at a radius of about 1.4 x 10^* cm (or 16" for HD 179821) from the star (Mamon, 
Glassgold & Huggins 1988). Thus, in agreement with observations, the models predict the 
CO persists to much greater distances from the star than does the OH. Once the CO is 
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photodissociated, the atomic carbon is rapidly photo-ionized by the ambient interstellar 
ultraviolet. As a consequence, little neutral carbon is expected to be found, consistent with 
the upper limit for the C I emission reported by Knapp et al. (2000). 

The observed circumstellar envelope around HD 179821 is consistent with that expected 
from theoretical models for the evolution of massive stars. In the computations by Schaller 
et al. (1992), stars with initially 25 Mq and 40 Mq have maximum mass loss rates as red 
hypergiants of 3 x 10~^ Mq yr~^ and 3 x 10""^ Mq yr~^, respectively. In the calculations 
by Chin & Stothcrs (1990), the mass loss rate for a star of 30 Mq at times reaches 4 x 10^^ 
Mq yr^^, the star shrinks to approximately ~10 Mq and the duration of intense mass loss 
phase is ~8 x 10^ yr. 

4.2. Deviations from Spherical Symmetry 

While the spherically symmetric model can reproduce much of the data, it fails to explain 
all the observations. Consider the apparent size of the CO shell as a function of gas velocity. 
For a thin shell or radius, Rq, expanding at uniform velocity, \4o, the projected distance 
from the center of the nebula, AR, at velocity offset from line center, AV, is (Olofsson et 
al. 2000): 

AR/Ro = (1 - [AV/V^ry/' (2) 

Figure 6 shows a comparison between the full width half maximum of the azimuthally- 
averaged CO intensity in the different velocity intervals. The model is a reasonable first 
approximation to explain the data, but there may be deviations by a factor of 20% in the 
outflow velocity from that predicted by this model. 

Jura & Werner (1999) found that HD 179821 does not he at the center of the circum- 
stellar ring, and a plausible interpretation of their data is that the outflow speed varies by 
as much as ±20% from its average value. Also, as shown in Figure 5, the profiles of the 
blue and red components of the circumstellar CO fine at the center of the nebula disagree 
with each other, and therefore the outfiow kinematics are not the same in the two directions. 
Also, neither the blue nor the red component of the line profile can be reproduced by a single 
outflow velocity. 

To interpret the anisotropy in the CO map and to learn about the deviations from 
spherical symmetry, wc model the excitation of the CO molecule. First, we argue that 
the gas kinetic temperature is likely to be above 20 K. Following Jura, Kahane & Omont 
(1988) and Kastner (1992), the gas temperature is determined by a balance between heating 
caused by the supersonic streaming of dust grains through the envelope and cooling resulting 
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both from radiation and the adiabatic expansion of the matter. In the outer portion of the 
circumstcllar shell, where adiabatic cooling is more important than radiative cooling and 
where MV^ ~ L^/c, the gas temperature, T^, is given by the expression: 

Tk ^ Wgrain ngrain/nn] Q^'"^ {L^/c) (47r/Cij)"^ (3) 

where {cr grain grain /n^n) denotes the grain number density times the average geometric grain 
cross section divided by the number density of hydrogen nuclei, Q is the ratio of the effective 
cross section for momentum transfer of a grain compared to its geometric cross section and 
is Boltzmann's constant. Plausible but uncertain values oi Q = 0.01 and {(J grain ^ grain I ^h) 
= 10-2^ cm2 (Jura et al. 1988) yield Tk = 220 K. 

Since it is likely that Tk > 20 K in the regions which produce the CO emission, then 
the inferred mass loss rate is relatively insensitive to the gas temperature and the unknown 
spatial variations in the grain properties. Instead, Tb in a particular velocity interval at any 
offset angle varies as the rate of coUisional excitation or as the square of the local density, so 
Tb scales as (M/Voo)'^ (see equation 4 below). Therefore, the observed factor of 2 variations 
in Tb at offset angle of 4" are consistent with Voo and M deviating in different directions by 
as much as a factor of 1.4 from maximum to minimum. According to Jura & Werner (1999), 
a possible explanation for these variations is the presence of strong magnetic fields within 
the atmosphere of the mass-losing star (Soker 1998), but another model is a highly dipolar 
global convective flow (Jacobs, Porter & Woodward 1999). Magnetic fields sufficiently strong 
to affect the mass outffow morphology have been reported around the red hypergiant VX 
Sgr (Trigilio, Umana & Cohen 1998). Current models for magnetized winds from red giants 
(Langer, Garcia-Segura & Mac Low 1999, Garcia-Segura & Lopez 2000) do not reproduce the 
asymmetry that we observe around HD 179821; the observed circumstellar envelope around 
this star may provide useful constraints on future calculations. 



5. COMPARISON OF THE NEBULA AROUND HD 179821 WITH 

KEPLER'S SNR 

We now consider the possible fate of the nebula around HD 179821 when the star 
explodes as a supernova, and we compare it to Kepler's remnant which is only 400 years 
old. It has been previously suggested that the anisotropy in Kepler's remnant results from 
a bow shock produced by a wind from the pre-supernova (Bandiera 1987). Evidence for a 
pre-supernova wind is given by the knots of optical nebulosity in the circumstcllar remnant 
with proper motion speeds near 100 km s^^ (Bandiera & van den Bergh 1991). (This gas 
moving at ~100 km s~^ may have been partly accelerated by the supernova explosion and 



- 12 - 



thus have been produced by a wind with a speed near 30 km s~^.) Here, instead of a model 
with a bow-shock, we suggest that the anisotropics in Kepler's supernova remnant can be 
understood if the star exploded into an intrinsically anisotropic circumstellar wind similar 
to that currently found around HD 179821. 

There is a qualitative similarity between the 11.7 /im and CO maps of the shell around 
HD 179821 with the 21 cm continuum VIA map (Reynoso & Goss 1999) and the ROSAT 
map (Hughes 1999) of soft X-rays {E < 2 keV) of Kepler's supernova remnant. Both objects 
show very approximately circular symmetry, but with one side being markedly brighter than 
the opposite. However, there is also a quantitative difference between the two structures. 
In HD 179821, the contrast between the integrated intensity of CO to the north and south 
is about a factor of 1.4. The contrast of the peak intensity between the northern and 
southern hemispheres in the ROSAT image of Kepler's supernova remnant is approximately 
a factor of 7. Even though the X-ray images of Kepler's supernova remnant display a much 
greater contrast in intensity than does the CO image of HD 179821, below we suggest that 
the nebula-hke that around HD 179821 may evolve into an object hke Kepler's supernova 
remnant. 

A full hydrodynamical calculation of a spherically symmetric supernova explosion into 

an anisotropic prc-supcrnova wind is beyond the scope of this paper. Here, we consider 
a simplified version of the interstellar bubble models presented by Weaver et al. (1978). 
Let the supernova begin as a spherically symmetric explosion of energy Eq and propagate 
without significant radiative losses into a detached envelope whose density, p{R), is described 
as 

p = M/iAnVooR^) (4) 

for R > Rdetached and p = for < Rdetached whcrc Rdetached dcuotcs the inner boundary of 
the detached circumstellar shell. With the very approximate assumption that most of the 
mass in the observed remnant is swept- up material, Msw, then: 

1/2 M,^ R^ ^ Skin Eo (5) 

where ekin denotes the fraction of the supernova's explosive energy that goes into kinetic 
energy (instead of heating the gas). We may take Ckin = 0.4 (Weaver et al. 1978). In the 
hmit that R > Rdetached, then: 

Ms^^ MR{t)/V^ (6) 
The solution to equations (5) and (6) is: 

R = {9ekinEoV^)y'{2M)-yH'/' (7) 

If Eo = 10^^ erg, t = 400 yr and i? = 8 x 10^^ cm (Hughes 1999), then equation (7) 
yields M/K^ = 5.5 x 10^^ g cm^^ If Ko = 32 km s-^ then M was 3 x 10"^ Mq yr-^ 
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These parameters for the presupernova wind are moderately similar to those inferred for the 
envelope around HD 179821, which therefore may be following an evolutionary path parallel 
to that experienced by the star that produced Kepler's supernova. 

We now estimate the angular variation of the X-ray emission that might occur when 
the current shell around HD 179821 is shocked by a symmetric supernova explosion. When 
viewed tangentially, the intensity of the X-ray emission, Ix, can be written as: 



where p{Rimpact) is the density at the impact parameter, Rimpact, and f{T) is a complicated 
function of the gas temperature. At the expansion speed of Kepler's supernova remnant of 
1500-2000 km (Blair, Long & Vancura 1991), the post-shock gas temperature is ~ 3 x 
lO'^ K, and the emissivity of the gas near 1 kcV is relatively insensitive to the gas temperature 
(Raymond & Smith 1977). We therefore ignore variations in /(T), and from equations (4) 
and (8) write: 



of the X-ray intensity is about a factor of 3, which is somewhat smaller than the observed 
factor of 7 contrast in the X-ray emission. 

Chin & Stothers (1990) found that a star of initially 30 Mq on the main sequence 
shrinks to 11.3 Mq through the red hypergiant phase and then, after 76,000 yr, explodes as 
a supernovae. If this model applies to HD 179821, then during the next 76,000 yr, the shell 
expelled at 32 km s~^ will expand to an inner radius of 7.7 x 10^^ cm which is approximately 
equal to the radius of Kepler's SNR of 8 x 10^^ cm (Hughes 1999). 



Above, we have considered the effect of a symmetric explosion impacting upon an asym- 
metric envelope and the resulting supernova remnant. However, it is also imaginable that 
just as the outer material is ejected anisotropically during the slow wind phase, the inner 
material might similarly be ejected anisotropically during the supernova explosion. Such 
an asymmetric explosion could explain why some pulsars have space velocities near 700 km 
s~^ (Cordes & Chernoff 1998). If the total momentum loss. p. is expressed as the sum of 
Legendre polynomials, and if this expansion is simplified to just the first two terms, then: 



Ix oc 

fi!^) P^Rimpact) Rimpact 



(8) 




6. THE AFTERMATH OF THE SUPERNOVA 



p = [Po + Picos{e)][l/{A7r)]f 



(10) 
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where d is measured relative to the presumed axis of symmetry and r is a unit vector in the 
radial direction. The total z-component of the momentum carried away by the ejecta, p^, is 

(p-z) sin{e)de = (1/3) Pi (11) 



The magnitude of the total momentum carried away by the ejecta, pt, is 

pt = (27r) / (p-f) sin{e) de = Po (12) 
J 

If the supernova explodes with energy, Eq, and the material moves outward with velocity, 
Vexpiode, then Pt = 2EQ/Vexpiode- H thcic is a neutron star of mass, M„s, left behind after the 
supernova explosion, then the magnitude of the velocity kick, Vkick, acquired by the neutron 
star can be inferred from the assumption that the momentum carried by the neutron star is 
balanced by Pz, the momentum carried by the remnant, so that: 

Vkick = (l/3)Pi/M„, = (2/3)(P,/Po)Eo/{M^sVe.,i,ae) (13) 

Above, we have argued that the anisotropic outflow from HD 179821 is consistent with 
variations of M and Voo by 20% from the mean and therefore we adopt Pi/ Pq = 0.2. With 
standard values of Eq = 10^^ erg, Mns = 1-4 Mq and Vexpiode = 10^ km s^"*^, then equation 
(12) yields Vkick — 500 km s"^. Given that these parameters may vary by as much as a factor 
of 2 from the assumed values, it is possible that kicks as large as 700 km s"^ may result from 
supernova explosions if they are as asymmetric as the outflow from HD 179821. 

It is possible that supernova explosions produce gamma-ray bursts. If so, then the 
late-time light curves of these systems are in part determined by the propagation of the 
explosion into the pre-existing circumstellar nebula (Esin & Blandford 2000). When HD 
179821 explodes as a supernova, perhaps in 8 x 10^ yr, the inner radius of its current dust 

shell may have expanded to 2.6 pc. Esin & Blandford (2000) propose that supernovae possess 
dust shells at a distance in the range 0.1-1 pc to explain the 20-30 day delay in the dust 
echoes from gamma-ray bursts. The detached shell around HD 179821 may be representative 
of the environments where dust echoes from gamma-ray bursts are produced. 



7. CONCLUSIONS 

We have obtained high resolution J = (1— ^^0) CO observations of the circumstellar shell 
around HD 179821, a star which will probably become a supernova in the next 10^ yr. We 
make the following points: 
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• To reproduce the CO map and spectra, it seems that until about 1600 years ago, the 
star was a red hypergiant with an average mass loss rate near 3 x 10~^ M© yr~^ at an 
average outflow speed of 32 km s~^. 

• The north-south asymmetry in the CO emission resembles that seen at 11.7 iim. The 
CO data suggest that M/Voo varies by as much as a factor of 1.4 from minimum to 
maximum. There appear to be deviations by as much as factor of 1.2 in the expansion 
velocity of the gas. 

• The asymmetry seen in the infrared and molecular emission around HD 179821 resem- 
bles that seen in the radio and X-ray emission from Kepler's supernova remnant. Even 
if a supernova explosion is spherically symmetric, it might propagate into an intrin- 
sically asymmetric circumstellar nebulae and thus even young remnants can appear 
asymmetric. 

• If HD 179821 explodes as an asymmetric supernova with the momentum loss varying 
by as much as 1.2 in different directions, consistent with an extrapolation of our ob- 
servations of its outer circumstellar envelope to the inner core, then the pulsar that is 
created may have a space motion approaching 700 km s"^. 

• The detached dust shell around HD 179821 may be representative of the environments 
which produce dust echoes 20-30 days after a gamma-ray burst. 

We have had useful conversations with A. Kusenko and S. Phinney as well as email corre- 
spondence with R. Chevalier and E. Jossehn. This work has been partly supported both at 
UCLA by NASA and at the Jet Propulsion Laboratory, California Institute of Technology, 
under an agreement with NASA. The Owens Valley Radio Observatory Millimeter Array is 
supported in part by the National Science Foundation grant AST 96-13717. 
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FIGURE CAPTIONS 

Fig. 1. Plot of the LSR Velocity vs. the Outflow Velocity derived from CO data with 
squares for the oxygen-rich high mass loss rate (M > 10~^ Mq yr~^) AGB stars within 1 
kpc of the Sun from Jura & Kleinmann (1989). The crosses for HD 179821 and the well 
known hypergiant IRC+ 10420 are also shown. The data are consistent with the hypothesis 
that a high outflow velocity is correlated with a high luminosity, and that both HD 179821 
and IRC-l-10420 are distant, high luminosity objects. 

Fig. 2. The circumstellar dust and gas shells around HD 179821. The upper panel shows 
the J = 1-0 ^^CO intensity superimposed upon the 11.7 /im map obtained by Jura & Werner 
(1999); the lower map shows the CO data by itself. The spatial resolution of each map (l'.'4x 
2'/0 for CO and 0'.'4x 0'.'4 for 11.7 /xm) is shown. The CO intensities represent the integrated 
line intensities jTdv over a velocity spread of 65 km s~^ centered at Vj^,. = 100 km s~^. 
The first contour is at 142.4 K km s^-*^; the contours levels are spaced by 71.2 K km s^^ 
(1.8ct). The star's position is denoted as a red cross. The color-coding in the lower panel is a 
redundant indicator of the intensity of the emission. In both panels the highest intensity is 
represented by red and the lowest by blue. The image of HD 179821 in the 11.7 iim image is 
noncircular because of an imperfect chop-nod motion of the Keck telescope during the time 
that the data were acquired. 

Fig. 3. Channel maps of the circumstellar J = 1-0 ^^CO emission from the shell around HD 
179821 shown with 5.35 km s~^ velocity spacing. The star's position is denoted as a cross. 
The first contour intensity in each map is at 80 mJy beam~^ (2.6 K brightness temperature); 
the contours levels are spaced by 80 mJy beam~^ (2.3cr). 

Fig. 4. Azimuthally-averaged intensity of the CO emission (mJy/beam) in the channel 
map at fine center (0 km s~^) and the channel map at -|-27 km s"^ velocity offset as a func- 
tion of angular offset (") from the central star. 

Fig. 5. The CO spectrum (mJy/beam vs. velocity) at the center of the circumstellar 
envelope. 

Fig. 6. Plot of the FWHM of the CO emission (") for different velocity channels com- 
pared to the expectation from equation (2). For the model fitting, we take, the half width 
radius = 3'.'95 and Ko = 32 km s^^ The error bars in the FWHM of the CO emission are 
taken as 0'.'2 or 10% of the beamsize. 



40 



IRC+ 10420 



HD 179821 



30 

s 



I — I 

U 

o 
> 



20 - 



o 

D 10 







J L 



1 



J I I L 



1 



J I I L 







50 100 

|LSR VELOCITYl (km/sec) 



150 




Offset RA (arcsec) 




Radial Distance (arcsec) 




I I I I I I I I I I I I I I I I I I I 

2 4 6 8 10 

CO FWHM (arcsec) 



